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Abstract 

This paper is devoted to strictly hyperbolic systems and equations with non-smooth coefficients. 
Below a certain level of smoothness, distributional solutions may fail to exist. We construct gen- 
eralised solutions in the Colombeau algebra of generalised functions. Extending earlier results on 
symmetric hyperbolic systems, we introduce generalised strict hyperbolicity, construct symmetris- 
ers, prove an appropriate Garding inequality and establish existence, uniqueness and regularity of 
generalised solutions. Under additional regularity assumptions on the coefficients, when a classical 
solution of the Cauchy problem (or of a transmission problem in the piecewise regular case) exists, 
the generalised solution is shown to be associated with the classical solution (or the piecewise classical 
solution satisfying the appropriate transmission conditions). 

AMS 2000 MSGS: 35S30, 46F30; 35B65. 
Keywords: hyperbolic equations and systems; algebras of generalised functions; non-smooth coeffi- 
cients; symmetriser. 



Introduction 

This paper is devoted to generalised solutions of linear strictly hyperbolic systems and higher order 
equations with non-smooth coefficients. In particular, the coefficients might be so singular as not to 
admit solutions in the sense of distributions. Such a situation may arise when the regularity of the 
coefficients is lower than Lipschitz and the data and driving terms are distributions. Formally, this leads 
to products of distributions. Holder continuous coefficients or coefficients with jump discontinuities in 
space arise in wave propagation in non-smooth media and are of importance, e.g., in seismology [9]. 

In order to stay in a framework that admits solving hyperbolic equations and systems with strong coef- 
ficient singularities, we employ the Colombeau theory of algebras of generalised functions. Solvability of 
hyperbolic first order systems of differential and pseudodifferential equations in Colombeau algebras (or 
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their duals) has been estabhshed in the symmetric hyperbohc case [21, 27, 29, 30, 31]. What has been 
open is the strictly hyperbolic case and the case of higher order equations. One of the principal goals of 
this paper is to close this gap. For this purpose, we introduce the notion of generalised strictly hyper- 
bolic systems, construct a generalised symmetriser and prove inequalities of Carding type for operators 
with Colombcau symbols. This allows us to establish existence, uniqueness and regularity of solutions in 
Colombeau algebras. 

A second task is to relate the generalised solutions to classical ones, when the latter exist. Indeed, in 
certain sufficiently regular cases, classical solutions to the Cauchy problem do exist. For example, if the 
coefficients are piecewise smooth with jumps across smooth hypersurfaces, one may phrase the problem 
classically as a transmission problem (sec, e.g., [1] and references therein). In this case, two classical 
solutions may be constructed on either side of the interface and joined by suitable jump conditions. 
Further, in various hyperbolic systems in divergence form (e.g., the linearised Euler system in acoustics), 
i°°-regularity of the coefficients and iJ^- regularity of the data suffice for the existence of a classical 
solution. In such cases, we set out to show that - as a rule - the generalised solution is associated 
with the classical solution (i.e., the representing nets of smooth functions have the classical solution as 
their distributional limit). In the piecewise smooth case, the nets representing the Colombeau solution 
thus already encode a transmission condition. In physical systems (e.g., acoustics), convergence to the 
piecewise classical solution satisfying the physically meaningful transmission conditions will be shown. 
Simple results of this type have already been obtained in some special cases [27, 29, 30]. We will complete 
the known results in this direction by studying the linearised Euler equations in acoustics and the wave 
equation with discontinuous coefficients. 

To set our work in perspective, we give a short survey of classical results on hyperbolic equations with 
non-smooth coefficients. Systems with discontinuous coefficients have already been considered in the 
1950ies and 1960ies, using jump conditions and energy methods [17, 23, 26]. Second order hyperbolic 
equations with Lipschitz coefficients that depend on time only were a subject of research starting in the 
1970ies and 1980ies, culminating in the result that the Log-Lipschitz property constitutes the lower bound 
on regularity for having existence and uniqueness of a classical solution [6, 7]. For coefficients depending 
both on time and space, see also [5]. For the general theory of pseudodiffcrcntial and paradifFerential 
operators with symbols of low regularity we refer to the monographs [20, 38] . Another road of investigation 
has been the case of transport equations with discontinuous coefficients and/or measures as initial data. 
Renormalised solutions have been introduced in [11], _BF- vector fields have been studied in [2, 8], measure 
theoretic concepts have been exploited in [4, 33] . An exploration into products of distributions is in [22] . 
For an in-depth survey on these and further results we refer to [3, 19]. In summary, this line of research 
has aimed at pushing the regularity to its lower bounds in special equations. That there are lower limits 
is shown, e.g., by an example of a first order scalar equation in divergence form with a piecewise constant 
coefficient that does not admit a solution in the sense of distributions, even for constant initial data 
[23]. In contrast to this, we place ourselves in a general framework that admits generalised sohitions for 
arbitrary (strictly) hyperbolic equations and systems with generalised function data. The relation to the 
previous literature is obtained by studying the limiting behaviour of the representing nets. For example, 
non-uniqueness of classical solutions is reflected in the dependence on the choice of rcgularisation. In 
ongoing research [10, 15, 16, 31], we are also developing an intrinsic characterisation of the properties of 
the solution in case no limit exists, exploiting the asymptotic behaviour of the representing nets. 

The plan of the paper is as follows: we begin by recalling some notions from the theory of Colombeau 
algebras as well as on Sobolev continuity of pseudodiflFerential operators. In Section 2, we consider linear 
hyperbolic systems with Colombeau symbols, introduce the corresponding notion of strict hyperbolic- 
ity and construct a symmetriser. The symmetriser is represented by a net of smooth functions whose 
asymptotic bounds have to be analysed so that the right mapping properties between Colombeau alge- 
bras are obtained. Next, we show that the symmetriser can be chosen positive definite in an appropriate 
generalised sense. To perform the symmetrisation up to a regularising error, the asymptotic scale of 
the rcgularisation parameter has to be chosen appropriately (slow scale to have existence of generalised 
solutions, logarithmic slow scale to have regularity, in addition). The necessity for these types of scales 
has been pointed out, e.g., in [14]. In Sections 3 and 4 we establish existence, uniqueness and regularity 
results for strictly hyperbolic systems and higher order equations, respectively, in the Colombeau setting. 
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Coming naturally with energy estimates, the results are formulated in terms of the Colombcau algebra 
based on Sobolev spaces. Section 5 addresses the limiting behaviour in the case of sufficiently regular 
coefficients. In the system of miiltidimcnsional acoustics, we show convergence to the classical sohition for 
L°°-cocfficients and H^-data. For the one-dimensional wave equation with piecewisc constant coefficients, 
we show convergence to the distributional solution of the corresponding transmission problem. We em- 
ploy Kato's perturbation resiilt for semigroups as well as Sobolev space techniques. Finally, the Appendix 
contains a proof of the required Carding inequalities for generalised pseudodifferential operators. 

1 Basic notions 

This section collects some preliminary notions concerning the different kinds of nets and of quotient spaces 
used in the paper. For the sake of brevity we mainly report definitions and basic properties, referring to 
[13, 14, 18, 34] for further details. 

1.1 Colombeau theory 

Nets of numbers 

A net (Me)^ in C*^°'^l is said to be strictly nonzero if there exist r > and r] G (0, 1] such that \us\ > e'" 
for all £ G (0, 7]]. 

For several regularity issues we will make use of the following concepts of slow scale net and logarithmic 
slow scale net. A net € M^'-^l is a slow scale net iff \oJe\ = 0{£~^) for all p > 0. We say that G M°'-^l 
is a logarithmic slow scale net if \uj^\ = 0(log^(l/e)) for all p > 0. Note that the suffices sc and Isc will 
stand for slow scale and logarithmic slow scale respectively. 

Generalised functions based on a locally convex topological vector space E 

The most common algebras of generalised functions of Colombcau type as well as the spaces of generalised 
symbols we deal with are introduced by means of the following general models. 

Let E he a locally convex topological vector space topologised through the family of seminorms {pi}iei- 
The elements of 

VielBNGN pi{ue) = 0{£-^) ass 0}, 
Vi e / 3(we)£S.s.n. Pi{ue) = 0{uis) as e 0}, 
3N€N\/i€l K(iie) = Ci(e-^)ase ^ 0}, 
Vi e J Vg e N pi{us) = 0{e'i) as£ ^ 0}, 

are called £^-moderate, £^-moderate of slow scale type, £^-regular and negligible, respectively. We define 
the space of generalised functions based on E as the factor space Qe '■= Me/J^e- 

The ring of complex generalised numbers, denoted by C, is obtained by taking E = C We remark that 
C is not a field since by Theorem 1.2.38 in [18] only the elements which are strictly nonzero (i.e. the 
elements which have a representative strictly nonzero) are invertible and vice versa. Note that all the 
representatives of u € C are strictly nonzero as soon as there exists one representative which is strictly 
nonzero. 

For any locally convex topological vector space E the space Qe has the structure of a C- module. The 
C-module Ge ■— M.^eI-^e of slow scale regular generalised functions and the C-module Q'^ :— M.'^ /Me 
of regular generalised functions are subrings oiQE^ which are characterised by more specific moderateness 
properties at the level of representatives. We use the notation u = [{ue)e\ for the class u of {ue)e in Qe- 
This is the usual way adopted in the paper to denote an equivalence class. 



Me := {{ue)e G : 
Me := {iu,% e : 



3 



Colombeau algebras and spaces of generalised symbols used in the paper 

The Colombeau algebra Q(fl) on an fl open subset of K" is the factor space Qe with E = C°°{fl). For 
the hyperbolic problems studied in this paper we will mainly choose E = C^{[—T,T] x M"), the space 
of infinitely differentiable functions all whose derivatives are bounded, and work in the corresponding 
Colombeau algebra Q^, ( [— T, T] x R") = Gc^([-t.t] xR") • It will be also convenient to work with Colombeau 
generalised functions based on Sobolev spaces. In particular, t/2,2(IR") := 5/f°°(R") and G2,2{{—T,T) x 
M") := Qh'>={{-t,t)xR'^)- 

With the expression generalised symbol we mean an element of the Colombeau space Qe, where E is 
a suitable space of symbols. For example E = S""(R2")^ e = C°°([-r, T], ^'"(M^")), etc. Specific 
moderateness properties can be required at the level of representatives. We recall that the slow scale 
regular generalised symbols {a = [{as)s] € Q^m(jg2n-j) are those generated by nets in M^gm(^2„y more 
precisely, 

(ae)e e Al5m(R2n) -^=> for all a, /3 e N" there exists a slow scale net such that Ictel^"^^ = 0{LJe), 

where |a|;"''' denote the seminorm sup(^_j')gR2,i |£j_|_^|<; {^)-^+M\d^d^a{x,^)\. Analogously one can define 
symbols which are logarithmic slow scale regular by taking representing nets with the following property: 

for all a,j3 ^ N" there exists a logarithmic slow scale net u>e such that |ae|^"j3 = 0{lo^). 

We use the notations [a^)^ G A4'|Si|-R2r,-| and a e (]Tj2n ) . In [13, 14] a complete symbolic calculus has 
been developed for slow scale regular generalised symbols. In the same vein, a complete symbolic calculus 
can be developed for logarithmic slow scale regular generalised symbols as well. 

In general, the net of representing smooth functions {ue)e of an element u G Q{Vl) does not converge. 
However, if it converges weakly to a distribution v £ 2?'(^)i then u is said to admit v as associated 
distribution. 

1.2 Sobolev-boundedness of generalised pseudodifferential operators 

We collect some basic facts concerning pseudodifferential operators which will be employed in the sequel. 
Our main reference is [25]. 

Sobolev spaces 

We recall that {D^)'^ denotes the pseudodifferential operator with symbol {^Y = (1 + |^|^)* and that the 
Sobolev space H^{W) is the set of all u € ^'(K") such that {D^Yu e L^{M"). H^R") is a Banach space 
with respect to the norm \\u\\s = \\{Dx)'^u\\i2(^^,iy When u G iJ''(R") and v € iJ~''(M"), the inequality 
\{u,v))\ < ||u||s||w||_s holds for the distributional action of v on u. 

Sobolev-boundedness theorem for pseudodifferential operators on K" 

Let ^"'(M^") be the space of Hormander symbols of order m satisfying global estimates in x and ^. The 
corresponding pseudodifferential operator a{x,D) defines a continuous map from /f*+™(M") to iJ*(M"). 
More precisely, see [25, Theorem 2.7], for each s there exist a constant Cg and an integer Ig such that 

(1.1) \\a{x,D)u\\s<Cs\a\^i^'>\\u\\s+m 

for all u G Note that the constants Cg and Ig only depend on s and the order of the symbol 

a. 



4 



Sobolev-boundedness of nets of pseudodifferential operators and of generalised pseudodif- 
ferential operators 

Let us now consider a net of symbols {ae)e G . It follows from the above that the inequality 

(1.1) holds uniformly with respect to s: in other words, 

(1-2) \\a,ix,D)u\\s<C,\a4f\\u\\s+m 

for all e £ (0,1] and u G _ff^+'"(]R"). In particular, the net {Cs\a^\l"^'')e of continuity constants will 
inherit the moderateness property of the nets of symbols {ae)e- Hence, a generalised pseudodifferential 
operator a{x,D) with symbol a G ^5m(R2n) is well-defined as a map 

a{x,D) : (7^e+™(R„) -> GH^iW")- 

This map is continuous with respect to the corresponding sharp topologies on the Colombeau spaces 
Qfjs+m(^n-^ and Qh^(&") (see [12] for the definition of sharp topology). 

Remark 1.1. Let ^2,2(K"') be the Colombeau space based on H°°(R"-). Prom (1.2) it immediately 
follows that any generalised pseudodifferential operator a{x,D) maps ^2,2 (IR") continuously into itself 

2 Symmetrisable systems 

This section is devoted to first order systems of the following type: 
(2.1) ^u = Ku + f, 

where K(t, x, D^) is an m x m-matrix of generalised pseudodifferential operators kij (i, x, D^) with symbol 
in the space Qc°°{[-t,t],s^{m?'^))- We use the notation K{t,x,^) for the corresponding matrix of symbols 
with entries kij{t,x,^), i,j = 1, ...,111. We say that a G Ge°° {[-t.t],S"^{m^")) is a classical (or polyhomoge- 
neous) generalised symbol if there exists a sequence of symbols am-j G Sc°°([-t,t],S'"-j(r2")), a sequence 
of representatives {am-j,e)e and a representative (oe)^ of a such that the following holds: 

- am-j,e is homogeneous of degree m — j in ^ for |^| > 1, i.e., 

for ]^| > 1 and r > 1; 

- for each r the net (ae - Z^J=o O'm-o,e)e G ■A4c~([-t,t],S"'-'-MR^"))- 

We refer to am G 0c°°([-t,t].s™(R2")) as the principal part of a G Qc°°{[-t,t],S"^('s?"')) or principal symbol 
of the operator a{x,D). Clearly it is uniquely defined modulo 0c°°([-t,t],S'"-i(r2"))- Note that when one 
works with more specific kinds of moderateness (i.e., slow scale, logarithmic slow scale) the same kind of 
moderateness is involved in the corresponding notion of asymptotic expansion. Hence the principal part 
of a slow scale (or logarithmic slow scale) regular classical symbol is slow scale (or logarithmic slow scale) 
regular as well. 

2.1 Symmetriser: definition and existence in the generafised strictly hyper- 
bolic case 

From now on we assume that K(t,x,^) has a classical generalised symbol, i.e., with entries in the space 
&C~([-T,T],si(R2")) and we denote its principal symbol by Ki{t,x,S,). We introduce the following notion 
of symmetriser and symmetrisable system. 
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Definition 2.1. A symmetriser for d/dt — K is a smooth one parameter family of operators R = 
R{t, X, Dx) such that 

(i) the entries of R are given by classical generalised symbols in Qc°°([-t,t].s<^(K^^))'j 

(ii) Ro{t,x,^) is a positive definite matrix for |^| > 1, in the sense that there exists a representative 
{Ro,s)s and a constant c > such that 

z'^Ro^e{t,x,(,)z > c\zf 

for all z e C", for all t e [-T, T], x eW, e € (0, 1] and for all \^\ > 1 

(Hi) RK + [RK)* is a matrix of pseudodifferential operators with symbol in Ge°°([-T,T],s°{R'^")) ■ 

If such a symmetriser exists, we say that d/dt — K is symmetrisable. 

Definition 2.2. The operator d/dt — K is called generalised strictly hyperbolic if there exists a repre- 
sentative {Ki^g)g of the principal symbol Ki{t,x,^) having m pure imaginary eigenvalues {i\j^e)e o/nd a 
strictly nonzero net (Ae)^ such that the bound from below 

(2.2) A,+i,e(t,a;,0 - XjAt^x,^,) > X,{0 

holds for all j = 1, m — 1, for all e £ {0,1], t G [— T, T], x € R" and for ^ away from zero. 

The expression ^ away from zero stands for |^| > r with < r < 1. 

After truncation at zero, the eigenvalues iXj in Definition 2.2 generate symbols with the same asymptotic 
properties (in e) as the matrix Ki{t,x,^). Indeed, using arguments as in the proof of Proposition 6.4 in 
[28] one gets the following result. 

Lemma 2.3. The nets {Xj.e) generates m, distinct and real classical symbols \'j in Qc'>^([-t,t],S'^(m?^)) 
such that \'j{t,x,£/) = \j{t,x,£/) for \£\ large enough. 

Proof. Wc work at the level of representatives. Let Ps{\;t,x,^) := (det(A/+i/i'i,£(t, a;, ^)) = ni<j<TO(A — 
\jAt,x,^)). The bound from below (2.2) yields 



dPe^ 
dX 



X=Xi 



m— 1 



for t G [— T, T], X e R" and |^| > r. Since each eigenvalue can be bounded by the norm of the matrix 
KiAt,x,£/) we have the estimate 

|A,-,e(i,a;,OI <Me(C>, 

valid for t G [—T, T] and [x, G M^", where (/ie)^ is a moderate net depending on the coefficients of Ki^^. 
Passing to the first order derivatives from Ps{X;t,x,£/) = ni<j<m(A — \jAt,x,S,)) we get, for |^| > r. 



An induction argument shows that each eigenvalue \jAt,x,^) satisfies moderate symbol estimates of 
order 1 for |^| > r. Prom the fact that Pe is homogeneous of order m in ^ for |^| > r we obtain that each 
net {\j,e)e is homogeneous of order 1 on the same ^-domain. Let now ij) G C°°(R"), ^/'(^) = for |^| < r 
and V(C) = 1 for ICI ^ Q^''; '^ith cur < 1. The net Xj^^{t,x,(;) = ^jj{^)Xj^s{t,x,^) generates a real classical 
symbol A^ in Gc=°{[-t.t],s^{r^^)) which coincides with Xj for |^| > ar. □ 

Remark 2.4. From the previous lemma it is clear that when Ki has slow scale entries then the eigenvalues 
are slow scale regular as well. Note that one can slightly generalise Definition 2.2 by allowing a strictly 
nonzero radius {re)^. Finally, in order to get slow scale nets of eigenvalues one combines the slow scale 
assumption on the entries of Ki with the condition that (re)e as well as (Xe)e is the inverse of a slow 
scale net. 
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Remark 2.5. In the sequel we will make use of the adjoint A(t, x, D^)* of a matrix A{t,x,Dx) of 
pseudodifferential operators. By definition A(t, x, D^)* fulfills the relation {A{t,x, Dx)u,v)L2(jgin^m — 
(u, X, Z)a;)*w)i2(Rn)m for all u,v G ^(M"). At the symbol level we have that, modulo lower order 

T 

terms, A(t, x, ^)* = A{t, x, ^) . All these statements are valid for generalised pseudodifferential operators 
as well. 

Proposition 2.6. Any generalised strictly hyperbolic operator d/dt — K is symmetrisable. 

Proof. We have to construct an operator matrix R fulfilling Definition 2.1. We begin by writing the 
projections Pj,^ onto the associated eigenspaces of iXj^e- For the moment we take ^ 7^ and we work 
with ^/|^|. Since from the hypothesis of generalised strict hyperbolicity we have 

Aj+i,e(t,a;,^/|C|) - Xj,e{t,x,^m) > X„ 

with a strictly nonzero net {X^)e, we can find a curve jj^dt, x, C/ICI) around the point iXj^e{t, x, C/ICI) such 
that all iXh,£{t,x,^/\^\) with h ^ j do not belong to the closed domain with boundary 'yj^s{t,x,^/\^\). 
Hence, 

Easy computations, as in [36, Proposition 2] lead to the following formula 

-zi^i,e(i,x,^/|e|)-Aft,e(i,x,C/|^|)/ 



The hypothesis (2.2) of strict hyperbolicity combined with the fact that Ki^^ is homogeneous of order 1 
in ^ implies that Pj^e is homogeneous of order in ^ (for ^ ^ 0). More precisely, from Lemma 2.3 and 
(2.4) we easily obtain that Pj^g{t,x,^) = Pj,e(t, .x, ^/|^|) is a matrix with entries in A4e°°{[-T,T].s°{M^"\a))- 
Truncating around (for instance taking a cut-off function identically for |^| < 1/2 and identically 
1 for 1^1 > 1) and keeping the same operator notation we obtain that the operators Pj{t,x,Dx) = 
[{Pj^t, X, Dx))s] are defined by matrices Pj{t, x, ^) of classical generalised symbols in ^c~([-t,t],s<'(m2"))- 
Let now 

m 

R{t, X, D^) := J2 Pj{t> a;, D^)*Pj{t, x, D^), 

where for the definition of Pj{t,x, D^)* we refer to Remark 2.5. The entries of the corresponding ma- 
trix R{t,x,^) are classical generalised symbols in Gc°°{[-t,t].s°{s.'^'^))- Moreover, the principal symbol 

Ro(t,x,S,) is given by X^^i Pi(i, a^, C) Pjit^x,^). Choosing the representatives Pj^^ as above we have 
that Ro,e{t,x,^) = Y,]Li Pj,ei't, X, ^)* Pj^e{t, X, ^) fulfills the following bound from below for any ^ e C" 
and for '|^| > 1: 



z^Ro,,{t,x,^)z = \PjAt,x,Oz\'' > \ 



Before proceeding we observe for ^ ^ we have X^^^ «Aj(f, x, ^)P,(t, cc, ^) = iiri(i,x,^), Pj{t,x,^) = 
Pj{t,x,^) and PjPh{t,x,^) = for j ^ h. It follows that modulo Qc^ {\-t,t\,s°(m?'^)) the symbol of RK 
is given by 

m 

i [t, X, ^)Pj{t,x,if Pj (i, X, 0- 
i=i 

Hence, once more modulo ^c°°([-t,t],so(r2")) the symbol of {RK)* is 

m 

-iJ2^j (*' ^' OPj{t,x,0'^Pj {t, x,0- 

i=i 

and then (RK) + (RK)* is a matrix of generalised pseudodifferential operators of order 0. This means 
that i? is a symmetriser for the operator d/dt — K. 

□ 
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2.2 An L^-positive definite symmetriser and tlie corresponding L^-energy es- 
timates 

Note that the symmetriser R constructed above is L^-self- adjoint, i.e., R = R* . Moreover a straight- 
forward apphcation of the Garding inequahty as stated in the Appendix yields a symmetriser S which 
is positive definite on L^(IR"). In the sequel (•,•) denotes the scalar product in L^. The corresponding 
norm is denoted by || • || or cquivalcntly || • ||o (in the Sobolcv space notation). 

Proposition 2.7. Let R be the symmetriser constructed in the proof of Proposition 2.6. Then there 
exists a symmetriser S such that 

(i) R = S (modulo a matrix of pseudodifferential operators with symbol in Qc°^{1-t,t],s-'^(r^")))^ 
(a) S is L"^ -positive definite, i.e. there exists a representative {3^)^ and a constant c > such that 
(2.5) {SeU,u)>c\\uf 
for all u G i2(M")'" and all e e (0, 1]. 

Proof. We begin by estimating Re i?£(f,a;,^). Since 

Rei?e(i,a;,0 = -Ro,e(t,a;,0 H 

we obtain from the properties of i?o,e and the symbol order of R-\^e + R*-i,e ^l^^t 

This bound from below holds for t G [-T,T], (x,^) e M^" and s G (0, 1] with the moderate net {iJ,e)e 
depending on the symbol estimates of J?_i,£ + R'^i ^. Hence for := 2m^/i£ and for |^| > r^ we get 

z'^ReR,{t,x,Oz > -^\A^- 
It follows that for c = l/(2m^) the inequality 

Reii£(t,a;,0 > c/ 

holds for 1^1 > r^. An application of the Garding inequality in Theorem A.4 entails 
Re(i?E(t, X, Dx)u, u) > c\\u\\q — ci,£||u||^i/2' ^ i^(K"), 

for some moderate net (ci,e)e. Let us now define 

Se{t,x,D,) = R,{t,x,D,)+ci,e{D,)~^I. 

It is clear that the corresponding generalised operator S is equal to R modulo a matrix of operators of 
order —1. Furthermore it is L^-positive definite. Indeed, since Ss{t,x,Dx)* = Se{t,x,Dx) we have that 
{SeU, u) is real and from the Garding inequality above 

{SeU,u) = {RsU,u) + Ci^e{{Dx)~^Iu,u) > c\\u\\l ~ Ci^e\\u\\^^i/2 + Ci^e{{Dx)~^ Iu,u) 

= c\\u\\l - ciM-1/2 + ci,,{{Dx)-^/^Iu, {D,)-^I^Iu) 

= 4Al - Cl,e||M||-l/2 + ClAW\W/2- 

Concluding 

{SeU,u)>c\\u\\l. 

□ 



8 



The symmetriser S yields an immediate energy estimate for the net of solutions (ue)e to the equation 
(2.6) —u, = KS,x,D^)u, + f„ 

with initial data Ue{0,-) = g^- More precisely, we assume that the operator d/dt — K is generalised strictly 
hyperbolic, we take {Ke)e with principal part (ii'i,e)e as in Definition 2.2, e C°°{[-T,T], L'^{W))'^ 
and g, e L^iR"")"". We have 

= {dtS,U„ Ue) + {S,K,Ue + SJ„ U,) + {S,U„ K,U, + /e) 

(2-7) < c,{dtS,)\\u,\\^ + c,{SeK, + K*S,)\\uef + 2c,{SM\\\\u4 

< cAdtS,)\\u,f + C,{S,K, + KlSMUe\? + C.(5e)(||u.||' + ll/ell') 

where the net (Ce)e depends on symbol seminorms of {Ss)ei {9tSe)e and {SeKs+K*Ss)e and the constant 
c appears in the bound from below (2.5). 

An application of Gronwall's lemma gives the estimate 

(2.8) \K{t,-)f < C-'(5eWe,We) < C et) , 

vahd for all s e (0, 1] and t e [-T, T] with 

d, = {Se{0,x,D^)g„g,) + Ce sup ||/e(i, Oll'^^- 

-T<t<T 

2.3 i^'-energy estimates for the solution net {u^)e 

The i^-positivc definite symmetriser S constructed in the previous subsection can be employed to ob- 
tain if'-energy estimates for the net {ue)e for any I e M. Before proceeding with the proof we re- 
call that i?' (and analogously (iJ')™) is a Hilbert space with the scalar product defined by {u,v)i = 
{{DxYlu, (DxYlv), where {DxY is the pscudodiffcrcntial operator with symbol (^)'. For the sake of 
simplicity in the sequel we use the short notation {D^Y for the matrix {DxYl- We now argue on 
UeJe G C°°{[-T,T],H'{W))"' and g^ e il'(M")'". The fact that is i^-positive definite yields 

\\u,\\f = \\{D,Yuef < C-\S,{D,Yue, (D^Yu,). 

In analogy with (2.7) we have 

j^{Se{D,Yue, {D^Y^,) = {dtSe{D,Yue, {D^Y^e) + {S,{D^YdtU„ {D^Y^e) + {S,{D^Yue, {D^YdtUe) 
= {dtSe{D,Yue, {D^Yue) + {S,{D^YKeU, + S,{D^Y fe, {D^Y^e) + {Se{D,Yue, {D,Y{KeU, + fe)) 

< C,{dtS,)\\Ue\\f + {S,{D,YKeUe, {D^Y^e) + 2Re {S,{D,Y fe, {D^Y^e) + {Se{D,Yue, {D^YK.U,) 

< C,{dtSe)\\Uefi + Ce(5e)(||We||? + ||/e||?) + {Se{D,YK,U,, {D^Y^e) + {Se{D,Yue, {D^YKeUe). 

Writing {S,{D,YKeUe, (D^Y^e) + {SeiD^Yue, {D^YKeUe) as 

{S,K,{D,Yue, (D^Yue) + {Se[{D,Y,Ke]Ue, {D^Yue) 

+ {S,{D^Yue,K,{D,Yu,) + {S,{D,Yue, [{D^Y,K,]u,) 
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we arrive at 
(2.9) 

^(s,{D,yu,,{D,yu,) 

at 

< c,{dtS,,SMu,\\f + \\f,\\f) + {{S,K, + K:S,){D,yu,ADcoyu,) + 2Re{SSD,y,^ 

<Cu{dtS„s,,s,K, + K:s„s,[{D,y,KMK\\f + \\fs\\h 

< Cuc-HSe{D,yu„ (D^yu,) + Ci,4Mf. 

Returning to the Cauchy problem (2.6), we obtain by applying the Gronwall lemma to the energy esti- 
mates (2.9) that 

(2.10) \\ue{t,-)\\f < c~\Se{D,yu„{D,yu,) < c-'di,,eMc~'Ci,et), 

valid for all s e (0, 1] and t e [-T, T] with 

di,e = {Se{0,x,D,){D,yg„{D,yg,) + Ci,e sup \\f,{t, ■)\\'iT. 

-T<t<T 

and {Ci^e)£ depending on the symbol seminorms of {Se)e, {dtSe)e, {SeKe + K*Se)e and {SeHD^y, Ks])e- 

2.4 Choice of scale for the generahsed strictly hyperbolic case 

From the estimates (2.8) and (2.10) it is clear that the nets {c~^Ce)e and {c~^Ci,e)e need to be of loga- 
rithmic type in order to generate a Colombeau solution from the net {ue)^. This logarithmic behaviour is 
entailed by the following assumptions on the principal part Ki{t, x, D^) and the corresponding eigenvalues 
A,-. 

Proposition 2.8. Let Ki{t,x, D^) have symbol entries which are logarithmic slow scale regular. Let the 
operator d/dt — K{t,x, Dx) be generalised strictly hyperbolic with eigenvalues Xj such that the bound from 
below (2.2) is true for some net {Xs)e, inverse of a logarithmic slow scale net. Let S be the symmetriser 
constructed in Proposition 2.1. Then, the entries of S and SK + K*S are given by 0-order logarithmic 
slow scale regular generalised symbols. 

Proof. Under the assumptions that Ki j;{t, x, ^) and {Xe)e are both logarithmic slow scale regular one 
easily sees that the eigenvalues (Aj^e)^, as well as the projectors Pj^^lt, x, (^), j = 1, ...,m, arc logarithmic 
slow scale regular. It follows that the symmetriser R{t, x, D^), as constructed in the proof of Proposition 
2.6, is given by a matrix of generalised symbols which arc logarithmic slow scale regular, i.e. symbols 
in 0c^([^TT] ^"(R^"))- Hence RK + (RK)* has symbols entries in ^?c°<>([-tt] ^"(R^")) well. Following 
the construction of S from R as in the proof of Proposition 2.7 and by applying Garding's inequality of 
Theorem A.5(m) we easily see that S is generated by 

where {ci^e)e is a logarithmic slow scale net. Hence, it has symbol entries in ^c°<»([-TTj s^iR^"-))' ^ ~ 
SK + K*S is of order and logarithmic slow scale regular and {S^u, u) > c\\u\\ for all u G L^(M")'". □ 

Corollary 2.9. Let Ki{t,x, Dx) have symbol entries which are logarithmic slow scale regular. Let the 
operator d I dt — K(t,x, Dx) be generalised strictly hyperbolic with eigenvalues Xj such that the bound from 
below (2.2) is true for some net {Xe)e, inverse of a logarithmic slow scale net. Then in the estimate (2.10) 
one has {Ci^e)e of logarithmic slow scale type. 

Proof. {Ci,e)e depends on the symbol seminorms of {Se)s, {dtSe)e, {SeKf. + K*Ss)e and {SeiiDxY , Ks\)s. 
From Proposition 2.8 it follows that (Cj,e)e is a logarithmic slow scale net. □ 
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3 The Cauchy problem for a first order generalised strictly hy- 
perbolic system 



Wc arc now able to prove the well-posedness in the Colombeau framework of the generaUsed strictly 
hyperbolic Cauchy problem 

d 

■^u = K{t,x,D^)u + f{t,x), 
w|t=o = g, 

where 



(hi) K{t, X, Dx) is an m x m-matrix of generalised pseudodifferential operators with classical symbol 

belonging to 0c~([-t,t],si(r2")); 

(h2) / has entries in a2,2([-T,T] x M"); 

(h3) the operator d/dt — K{t, x, Dx) is generalised strictly hyperbolic with eigenvalues Aj such that the 
bound from below (2.2) is true for some net (Ae)^, inverse of a logarithmic slow scale net; 

(li4) g e g2,2(K")". 

Theorem 3.1. Under the hypotheses {hi) — (M) the Cauchy problem (3.1) has a unique solution u e 
Q2,2{(-T,T) X M")™. /// as well as the initial data g are G^^-regular then u e G^^O^'')"'. 



Proof. Working at the level of representatives we obtain a net (u^)^ E H°°{{—T,T) x M")™ of classical 
solutions to the Cauchy problem (3.1) (see [25, Chapter 7, Theorem 3.2] and [35, Chapter 5]). The 
Sobolev energy estimates (2.10) obtained by means of the symmetriser S in the previous section provide 
the desired Colombeau moderateness (for the i-derivatives one argues by induction directly on 9tUe = 
K^{t, X, Dj:)u^ + fs{t, xj). Finally, the net (ue)e generates a Colombeau solution u G G2,2{{—T, T) x M")™. 
The uniqueness can be easily checked via (2.10) in terms of negligibility estimates. Prom the energy 
estimates it follows that u e ^2^2 (^")™ if / ^^'i 9 ^''^^ 0^2'''^^^'^^^ ■ ^ 



4 Higher order hyperbolic equations with generalised coeffi- 
cients: reduction to a generalised strictly hyperbolic system 

The aim of this section is to study the Cauchy problem for hyperbolic equations of the following type: 

(4.1) Lu = f, -g^u\t=o = gj+i, i = 0, ...,m-l, 

where 

(4-2) L= — -J2Am-j{t,x,D,)g- 

and each Am~j is a differential operator of order m — j with Colombeau coefficients (the Colombeau 
space type will be specified in the sequel). 

We begin by performing a reduction to a first order system as in [37]. Let {D^) be the pseudodifferential 
operator with symbol (^). The transformation 
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with j = 1, ...,m, makes the Cauchy problem (4.1) equivalent to the following system 



(4.3) 



where 



d_ 
dt 



( «i \ 



( ip^) ... \ / "1 \ 

(p^) ... 

m 

\ h\ b2 bm J \ Um. J 



+ 



[f J 



with initial condition 
(4.4) 



bj = Am-j+i{t,x,Dx){Dx) 



Uj\t=o = {D^)"" ^Qj, j = l,...,m. 



The generalised strict hyperbolicity of the operator L in (4.2) can be equivalently stated in terms of the 
roots T{t,x,£,) of 

m— 1 

P{T,t,x,0 := {irr - E Am-j{t,x,0{iry, 
where Am-j[t,x,£,) is the principal part of Am-j{t,x,^), or in terms of the eigenvalues of the matrix 

B{t,x,0 ■■= 

where bj{t,x,0 = Am-j+i{t,x,0{0'~'^- 

Definition 4.1. The operator L in (4.2) is called generalised strictly hyperbolic if there exists a choice of 
representatives of the coefficients in Q{[—T,T] x R") such that the corresponding polynomial Ps{T,t,x,^) 
has m distinct real valued roots {Tj^e)s such that 



( 







. 


. 


\ 








(0 • 


. 












• (i) 




V 


bi 


6~2 




• bm 


J 



(4.5) 



Tj+i,^{t,X,^) - Tj^s{t,X,^) > TeiO 



holds for some strictly nonzero net {ts)^, for all j = 1, ...,m — 1, for all e G (0, 1], t £ [—T,T], a; e M" 
and for ^ away from zero. 



For the same choice of the coefficient representatives we have that the purely imaginary nets (iTj^e)^, 
j = 1, m are eigenvalues of the matrix defined above and that Definition 4.1 coincides with Definition 
2.2. 

The well-posedness result of Section 3 for the Cauchy problem (4.3)- (4.4) yields the following statement. 
Note that for the Sylvester matrix in (4.3) we use the set-up of Theorem 3.1. 

Theorem 4.2. Let L he a differential operator of order m as in (4.2) with coefficients in Qb{[—T, T] x R") 

such that its symbol belongs to ^c=°([-tt] S"»(r2»i))- Let f G Q2,2{{—T,T) x E") and gj+i G ^2.2(R") for 
j = 0, ...,m— 1. Let L be generalised strictly hyperbolic with roots tj satisfying the bound from below 
(4.5) for some net (Te)e, inverse of a logarithmic slow scale net. Then there exists a unique solution 
u G a2,2((-T,T) X M") to the Cauchy problem (4.1). FuHhermore, if f e g^2{{-T,T) x M") and the 
initial data g^+i belong to 0^21^'") then u G QT^2{{-T-,T) x R"). 

Remark 4.3. The theorem holds, in particular, for logarithmic slow scale regular coefficients that are 
constant for large x as well as for logarithmic slow scale regular coefficients depending on t only, or 
belonging to a2,2((-T,T) x M"). 
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5 Applications and limiting solutions 



In the previous sections we have obtained gonerahscd sohitions to the Cauchy problem for hypcrbohc 
equations and systems with generahsed coefficients. As mentioned in the Introduction, classical solutions 
to the Cauchy problem may exist in certain sufficiently regular cases. It is the purpose of this section to 
show how the generalised solution in the Colombcau algebra relates to such a classical or picccwise classical 
solution, if it exists. We will show exemplarily (for the multidimensional system of linear acoustics and 
the wave equation) that the generalised solution is associated with the classical solution or the piecewise 
classical solution, if the appropriate transmission conditions arc imposed. I fact, in physical systems (e.g., 
in acoustics), the Colombeau solution encodes already the physically meaningful transmission condition 
and selects the corresponding piecewise classical solution in the limit. A similar result for a boundary 
value problem in one-dimensional acoustics has been obtained in [29]. 

As one of the central examples of interest in applications we consider the acoustics system. It describes 
small perturbations from a state of rest of the velocity v{x,t) £ M", the pressure p{x,t) € M and the 
density p{x,t) G IR of a substance in space {x G M") that moves in a medium of given density po{x) and 
sound speed co{x). It can be derived from the Euler equations of isentropic gas dynamics, written with 
material time derivative [39] . We will be particularly interested in the case where po and cq have a jump 
across some hypersurface, signifying a change of material properties at a separating interface. 

The acoustics system is 

d 

—p + Poix)divv = 
d 

(5.1) po{x)—v + gvadp = 

P = cl{x)p 

The divergence and gradient operators are understood with respect to the variables a; G IR". We are 
assuming that po and cq do not depend on t and are bounded away from zero. Thus we can cast the 
system in the form 



(5.2) 



—p + cl {x)po {x) div = 

—V + —— gradp = 
at po{x) 



We adjoin initial conditions 

(5.3) p{x,0) =pq{x), v{x,0) = vo{x). 

For sufficiently smooth solutions, the system (5.2) with initial conditions (5.3) is equivalent to the wave 
equation 

& 



(5.4) -^p - cl{x)po{x) div (^^^ gradpj = 



with initial conditions 



d 

(5.5) p{x, 0) = po{x), g^P^^' ^) = -Co(a;)po(a;) div vo{x). 

Indeed, differentiating (5.2) and subtracting the second from the first line immediately gives (5.4). For 
the other direction, one sets 

d 1 /■* 1 

— = ;r^gradp, v{x,t) = vo{x) + — --^ giad p{x,s)ds 

dt po{x) Jo Po{x) 

and uses the second initial condition (5.5) to deduce the first line in (5.2) from (5.4). 
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We will sec that the equivalence of the acoustics system and the wave equation pertains for solutions 
in the Colombeau algebra, hence we may freely switch from one formulation to the other. We will first 
consider the acoustics system (5.2) in any space dimension and then derive some more specific results for 
the wave equation (5.4) in one space dimension. 



5.1 The acoustics system in any space dimension 



Writing the system (5.2) in the form 



d 



—u = K{t,x,Dx)u, u\t=o=uo 
for the lumped vector u= {p,vi, . . . , Vn), we find that the symbol matrix is 



K{t,x,0 



-ipo{x)cl{x)^i 
a 



Po{x) 







-ipo{x)cl{x)^„ 








Its eigenvalues are A = ±ico (a;) a/^i + h and A = 0, the latter with multiplicity n — 1. Thus the 

system is not strictly hyperbolic for n > 2. For n < 2, Proposition 2.6 yields symmctrisability in the 
Colombeau setting. For arbitrary n > 1, it has been shown in [27] that the system can be symmetrized 
by a change of dependent variables. Nevertheless, in order to remain in the framework of the previous 
sections, we will turn to the wave equation formulation (5.4) for further analysis. Writing (5.4) as 



P - cl{x)Ap - cl{x)po{x) grad 



' gradp = 



we see that the principal symbol P{T,t,x,^) has the roots T(f, x, ^) = ±co(a;)|^|, which are real and 
distinct since co was assumed to be bounded away from zero. 

The results of Section 4 can now be applied to the wave equation (5.4) and the acoustics system (5.2). 

Theorem 5.1. LetcQ,po € Ghi^") be logarithmic slow scale regular generalised functions, hounded from 
below by the inverse of a logarithmic slow scale net. Let po,vo € G2,2(^") and T > 0. Then the initial 
value problem (5.4), (5.5) has a unique solution p G Q2,2{{—T,T) x R"), and the initial value problem 
(5.2), (5.3) has a unique solution {p,v) € a2,2((-T,r) x M")"+i. 



Proof. The existence of the solution p to (5.4), (5.5) follows immediately from Theorem 4.2. Under the 

assumptions on cq, po, it is clear that the acoustics system (5.2), (5.3) is equivalent to the wave equation 
(5.4), (5.5) in a2,2((-T, T) x R"). Thus the second assertion holds as well. □ 

Remark 5.2. To make the connection with possible classical solutions, consider coefficients co,po ^ 
L°°(IR") which are essentially bounded away from zero by some positive real constant. Take a compactly 
supported net {<fe)e of nonnegative functions which is logarithmic slow scale regular and converges to 
the Dirac measure. Then coe = cq * (fs and po,£ = Po* representatives of generalised functions 

Co, Po € Q{W^) that satisfy the hypotheses of Theorem 5.1. 



In order to prove the convergence of the representatives of the generalised solution to a classical solution 
(when it exists), we need to make a digression into the classical theory. In order not to overload the 
notation, we will use the letters v, p etc. in system (5.2) temporarily for classical solutions. We assume 
that Co and po are measurable functions satisfying bounds from below and above: 

(5.6) < pi < po{x) < P2, < Ci < Co{x) < C2. 
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We are going to show that the operator defined by (5.2) generates a semigroup of type (M, 0) where M 
depends only on the four bounding constants in (5.6). As usual, we let 

Fdi^(R") = e (l2(M«))" : divv e L'^iW')}. 

It is the domain of the divergence operator, which is a closed operator on (L^(R"))". We equip the 
Hilbert space 

X = L'^{W) X (l2(M"))" 
with the equivalent inner product (•!•) arising from the weighted norm 

llb,^)ll^o.o= / b(^)l %2J, M dx+ [ Hx)\'po{x)dx 
The unbounded operator A on X is defined by 

D{A) = H^W) X ifdiv(M"), A{p, v) = {-cIpo divv, -p^' gradp). 
In this setting, problem (5.2) reads ^(p, v) = A{p, v). The following assertion is the heart of the matter. 
Proposition 5.3. The operator A generates a contraction semigroup {S{t) : t >0} on X. 

Proof. The proof is standard, so wc just sketch the ingredients. First one shows by rcgularisation and 
cut-off that (C^(R")) is dense in iJdiv(R"), equipped with the graph norm of the divergence operator. 
In particular, A is densely defined. It is also quite obvious that A is closed. Next, 

(5.7) / {p{x) divt;(x) + gradp(x) • v{x))dx = 

for {p,v) G D{A). This follows from the fact that the left-hand side in (5.7) equals dW{p{x)v{x))dx, 
which vanishes for {p,v) e (C^(R"))""^\ 

The goal of the proof is to show that the operator A is maximal dissipative on X. As a first step, we 
observe that A is conservative, that is, {Au\u) = for all u = {p, v) e D{A). Indeed, 

{Au\u) = — I c^{x)pq{x) d\Y v{x) 2/ N — -i-^dx— f — ^-^ gi&dp{x) • v{x) po{x) dx = 

by (5.7). In particular, A is dissipative, i.e. (^u|w) < for all u = {p,v) G D{A). The fact that A is 
skew-adjoint, that is, D{A*) = D{A) and A* = —A follows from standard arguments and will be omitted. 
Thus A* is conservative as well, hence dissipative. Recall that an operator B is dissipative if and only if 

II (A - B)u\\ > X\\u\\ for all u e D{B) and A > 0. 

This implies that (A — B) is injective and has a closed range for all A > 0. We apply this to A and A* 
and infer that (A — A) has a closed range and (A — A*) is injective. The latter implies that the range of 
(A — A) is dense. Combining these observations, wc conclude that the range of (A — A) equals X for all 
A > 0. This is well-known to be equivalent to the property that A is maximal dissipative. 

Further, {X — A) : D{A) X is bijective and we have the resolvent estimate 

||(A-yl)-^|| < ^ for all A >0. 
A 

The Hille-Yosida-Lumer-Phillips theorem (cf. e.g. [32]) implies that A generates a semigroup of contrac- 
tions on X. □ 
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In particular, this means that the operator norm ||>S'(i)||copo ^ 1 for all i > with respect to the 
weighted norm on X. We wish to compute the operator norm ||S'(f)|| with respect to the usual norm on 
i2(]R") X (L2(Mn))". It is clear that 

CoPo 

for some constants Ci, C2 depending only on Ci, C2, pi, P2- Prom there it follows that 



(5.8) ||5WII<g||5WI|copo<§=M 



for all i > 0. In particular, {S{t) : i > 0} is a semigroup of type (M, 0) on L'^{W) x (^^(M"))". 

Corollary 5.4. Let% G ^^(M"), tJo e i?div(K") and co,po as in (5.6). Then problem (5.2) has a unique 
solution ip.v) m C^{[0,oo) : L2(M") x (L2(M")») such that {p{t),v{t)) belongs to H\W) x i?div(K") for 
all t >0, namely (^{t),v{t)) = S{t){pQ,Vo)- 

Proof. This is standard semigroup theory [32]. □ 

We are now in tlic position to relate the generalised solution with the classical solution, when it exists. 

Proposition 5.5. Assume that co,Po ^ L°°(M") satisfy the bounds (5.6). Let co,po G ^b(K") be 
corresponding generalised functions as defined in Remark 5.2. Let Pq G i?^(]R"), vq S iJdiv(K") and 
PQtVq G 02.2 (IR") be corresponding generalised functions with representatives {po^,vo^) obtained by con- 
volution with a delta net. Let T > 0. Then the unique solution {p,v) G G2,2{{—T,T) x ]R.")"+^ to (5.2), 
(5.3) given in Theorem 5.1 is associated with the classical solution {p{t),v{t)) given in Corollary 5.4, that 
is, its representatives converge to the classical solution as £ — >■ 0. 

Proof. As can be seen from Remark 5.2, the representatives coe, pos of the generalised functions co,po 
satisfy bounds of the type (5.6) with ci, C2, pi, P2 independent of e, due to the nonnegativity of the 
moUifier. Taking cqe, poe as coefficients in (5.2), Proposition 5.3 guarantees that the corresponding 
operator generates of a semigroup {Se{t) : t > 0}. Due to the classical and the generalised uniqueness 
result, the solution (pe(-, i), Ug(-, t)) = S'e (t) (poe , '^'Oe ) is a representative of the generalised solution {p,v). 
The estimate (5.8) shows that the semigroups {Se{t) : t > 0} are of type (M, 0) on the common Hilbert 
space L'^{W) x (l,2(R"))", independently of s. 

We now invoke Kato's perturbation result [24] which says that under these conditions, Se{t){po,vo) 

converges to S{t){po, vq) imiformly for t in compact time intervals, for every (po, vq) G I/^(1R") x (L^(M")) , 
provided the resolvents (A — A^)^^ converge strongly to (A — A)~^ for some A > 0. 

The resolvent convergence can be seen as follows. First take u = {p,v) G D{A) = D{As) G H^{W^) x 
iJ(jiv(K"). Then we have for the i^-norm 

\\AsU-Auf= \{cl^pos -lpo)divvfdx + \{pq^ - p^^) gv&dp\'^dx. 

Due to the bounds (5.6) and a wellknown property of Priedrichs mollifiers, the difference of the coefficients 

in the integrals converges to zero at every Lebesgue point, hence almost everywhere. By Lcbesgue's 
theorem, \\AsU - AuW^ ^ as e ^ 0. Now let u = (p,w) G L'^{W) x {L'^{R"))" be arbitrary. Then 

||(A - A,)-^u - (A - A)-'u\\ = ||(A - A,)-'{u - (A - A,)(A - A)-^u) \\. 

Since (A — A)~^u G D{A), the second factor converges to zero by the first step. Since {Ss{t) : t > 0} is 
a semigroup of type (M, 0), the resolvent norm ||(A — Ae)^^|| is bounded by M/A, uniformly in e. This 
entails that ||(A — As)~^u — (A — yl)~^M|| converges to zero, as desired. 

Finally, the fact that Ss{t){poe,voe) converges to S{t){po,vo) as well follows readily from Kato's pertur- 
bation result for fixed arguments and the fact that ||S'£(f)|| < M for alH > 0. □ 
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Remark 5.6. The classical solution {p,v) obtained in Corollary 5.4 belongs to i/^(R") x H^ivi^") for 
all times t > 0. li in addition it is piecewise continuously differentiable, then one can show that it satisfies 
the classical transmission conditions across the hypersurface of discontimiity. More precisely, let F be a 
smooth hypersurface that splits M" into two parts 17- , r2_|_. Let be the unit normal to F. One can 
show the following: If {p,v) belongs to C^(ri_) © C^{il+), then p and the normal component v ■ v are 
continuous across F. The second equation in (5.2) then implies that also gradp is continuous, which 
is the relevant transmission condition in the wave equation formulation (5.4). 

5.2 The one-dimensional wave equation with discontinuous coefficients 

We turn to the wave equation in one space dimension, which - slightly more general as in the previous 
subsection - we take of the form 

(5.9) a{x,t)-^w--[b{x,t)-w 
with initial conditions 

d 

(5.10) w{x,0) = wo{x), —w{x,0)=wi{x). 

The functions a and b are assumed to be strictly positive. We will treat the two cases where a and b 
depend on a; or on t only and are piecewise constant. We begin with the case 

(5.11) aix,t).-aix) = [l- ^ J K^, t) . 6(.) = { ^ < 
where the constants a_, a+, 6_, 6+ are strictly positive. 

Remark 5.7. If the initial function wq is continuously differentiable and wi is continuous, and a, & 
are piecewise constant as in (5.11), then there is a unique distribution w G I?'(R^) with the following 
properties: 

- IiJ is a distributional solution to (5.9) on the open halfplanes {x < 0} and {x > 0}; 

- w gC(IRx [0, oo))) and wJ(-,0) = uJq; 

- ^weC(Rx [0,oo))) and ^w(-,0) =wi; 



the function x b{x)-§-'w{x,t) is continuous for almost all t>0. 



This is easily seen by transforming (5.9) to a first order system on either side of x = and joining the 
two parts of the solution by the transmission conditions that w (and hence ■§fW) as well as b-^w should 
be continuous across a; = 0. We call W the classical connected solution of the transmission problem for 
the wave equation (5.9). 



As in the previous subsection, we produce logarithmic slow scale regular generalised functions a and 
b by means of representatives = a* ip^ with a mollifier as in Remark 5.2, and similarly for b^. By 
Theorem 4.2, problem (5.9) with initial data wq,wi £ ^2,2 (IR) has a unique solution w in ^2,2 ((—?", T) xR). 

Specifically, we take initial data mJq and wi as in Remark 5.7, which in addition - for the sake of the 
argument - have compact support. With any moderate mollifier converging to the Dirac measure, 
wo,£ = wJo * V'e defines an element of 52,2(R), and similarly with wig. Let w G Q2,2{{—T,T) x R) be the 
corresponding generalised solution. It vanishes if x is outside some compact set K. 

Proposition 5.8. In the situation described above, the generalised solution w to (5.9) is associated with 
the classical connected solution W. 
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Proof. Let (we)e be a representative of the generalised solution that vanishes for x outside K. This and 
the smoothness implies that belongs to C°°([0,T) : H°°{R)). Thus we can use energy estimates. 
Multiplying the wave equation awgu — {^ex)x = by w^t and integrating we get the conservation law 

\ d f 

-- {a\wet\'' +b\w,^\^)dx = fortG[0,T). 
z at 7k 

Differentiating the wave equations once with respect to time, multiplying by Wstt and integrating we get 




^)dx = iorte[0,T). 



The first conservation law implies that the family (w^)^ is bounded in C([0,T) : H^{R)). This implies 
that {we{-,t))e is relatively compact in C(M) for every t. Using the estimate 

\Ws{y,t) -We{x,t)\ < \/\y-x\ (^J \-^'l"e{y,t)\dx 

we get the equicontinuity of {w^{-,t))^. By Ascoli's theorem, (w^)^ is relatively compact in C([0,T) : 
C(]R)) = C(M X [0,T)). In the same way, the second conservation law shows that {'Wet)e is relatively 
compact in C(M x [0,T)). Further, {wett)s and {wsx)s are bounded in L^(1R x (0,T)). Thus for a suitable 
subsequence, we have that there are W,Wi gC{Rx [0,T)) and W2, W3 e L'^{R x (0,T)) such that 

We W, Wet -> Wi Strongly in C(M x [0,T)) and 
Wett W2, Wex W3 weakly in L^{R x (0,T)). 



It follows that Wi — Wt, W2 = Wtt and W3 = Wx- The wave equation now yields that {heWex)x converges 
weakly to some W4 G I/^(M x (0, T)). Due to the support properties, the convolution V{-,t) = Wi{-, t)*H 
exists, where H denotes the Heaviside function. Then V{-, t) is continuous for almost all t, and h^Wex V 
weakly in L'^{R x (0,T)). It follows that 

aWtt - K = 0. 

On the other hand, be{x) = 5+ at every given x > and he{x) = 6_ at every given a; < for sufficiently 
small £. Thus 

b-Wx, X < 0, 1/ / +^ _ / b-Wxx, a; < 0, 



\b+wZ x>Q, ^"^''^^ \b+Wxx, x>Q 

for almost all t. It follows that W solves the wave equation on either side of x = 0, W and Wt are 

continuous, take the given initial values uJq and wi, and bWx{-,t) is continuous for almost all t. By 
uniqueness, W is equal to the classical connected solution wJ, and the whole net We converges to it. □ 



As a final case, we consider the wave equation (5.9) with coefficients that have a jump in time, say at 
t = \. Thus we let 

(5.12) aix,t)^m = { ; 1 1; b{x,t)^m = [\- \l\ 

It is obvious that, given initial data 2?'(R), there is a vmiquc distribution w e C^([0, 1] : 2?'(R)) © 
C^([l,oo) : X''(K)) that satisfies the equation on either side of t ~ 1. As in the previous case, energy 
estimates and an application of a suitable negative power of the operator (D), if required, show again 
that the generalised solution is associated with the classical distributional solution. 
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A Appendix: Garding's inequality for a matrix of generalised 
pseudodifferential operators 

In the sequel we make use of the following notations concerning matrices of symbols: 

- Re^(x,0 = ^^^^^^^^^, 

- A{x, > if and only if z'^A{x, > for all zeC^ with z^O, 

- A{x,^) > B{x,^) if and only if {A - B){x,^) > 0. 

The proof of the Garding's inequality as stated in [25, Theorem 4.4] requires the following preliminary 
results on the Priedrichs part of a pseudodifferential operator. For more details we refer to [25, Chapter 
3, Section 4]. 

Definition A.l. Let q{a) G C^(|cr| < 1) be an even function with J q{a')'^ da = 1. Let 

Let p e 5'"(M^"). The Friedrichs part pF{i,x' ,^') of the symbol p{x , ^) is defined by 

Pf{^,x',0 = I F{^,0p{x',0F{^',0dC. 

Note that if p G S^CM^") then Pf{^,x',^') G i.e., 

\dp^!di'pp{^,x',a\ < cio—^i-'ir)"^'"'' 

for all a,a',/3' e N". 

Theorem A. 2. Let {pe)£ G Msm^M^ny Then 

(i) the corresponding net {pF,e)e belongs to A^^m.o ^^^a^y 

(ii) the operator pF,e{Dx,x' ,Dxi) can be written as a pseudodifferential operator with symbol {(JF,e)e € 

(Hi) [pe - (TF,e)e € Msm-i(^2ny 

The proof of Theorem A. 2 is obtained by combining Theorem 4.2 in [25, Chapter 3] with Theorems 2.5 
and 3.1 in [25, Chapter 2]. A careful inspection of the proofs shows that the expected moderateness 
estimates hold. Prom Theorem 4.3 in [25] we easily obtain the following statement. 

Theorem A. 3. Let Pe = Pe(x,Dx) be an I x I matrix of pseudodifferential operators with {Pe)e G 
A^gm(][t2n). Let PF,e be the Friedrichs part of P^. Then the operator Pf,£ — Pe is given by a matrix with 
symbol entries in Msm.-i(s_2ny In addition, 

(i) if Pe{x,^) is Hermitian, then PF,e is self-adjoint, i.e., 

{Pf,£U, v) = {U, PF,eV) 

for all u,v with components in ^(M"); 
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(ii) if Pe{x,£,) > then the operator Pp^e is positive, i.e., 

{PF,eU, U)>Q 

for all u with com,ponents in o5^(R"); 
(Hi) if Pe{x,^) is skew-Hermitian (P* = —Ps) then 

{PF,eU,v) = -{u,Pf,sV) 

for all u, V as above. 

We are now ready to prove the Garding inequality. For the sake of simpUcity we only consider 0-order 

generalised pscudodiffcrcntial operators. The norms used in the following theorem are the Sobolev norms 
of order and —1/2 respectively and (•,•) denotes the Hilbert product in L^(M")'" (or equivalently 



Theorem A. 4. Let A(x, D^) be anmxm matrix of generalised pseudodifferential operators with symbol 
in Qso(^2ny If there exist a representative {As)s of A{x,^) and strictly nonzero nets (ce)^ and (re)^ such 

that 

(A.l) ReAe{x,0>CeI 

for 1^1 > then there exists a moderate net (ci^e)^ such that 

(A.2) Re{As{x,Dx)u,u) > Cs\\u\\l - ci^s\\u\\ -1/2 

for all vectors u with components in and all e G (0, 1]. 



Proof Let Ai,s(a;,0 = Ae{x,C) - cj. Then 

Re{Ai,eix, Dx)u,u) = Rc{ A^ix, Dx)u,u) - Ci,\\u\\l. 

We define the operators Bs{x,Dx) and Ce{x,Dx) via the symbols ReAs{x,^) — c^I and {As{x,^) — 
As{x,$,)*)/2, respectively. Hence Ai^s{x,Dx) = Bg{x,Dx) + Ce(x, D^)- Making use of the hypothesis 
(A.l) and multiplying the entries of _Bg by a suitable cut-off net {il^e{0)s = {fp{£,/Te))e we obtain that 
there exist matrices £'£(3:,^) and Re{x,£,) with entries in A^50(r2„) and Ais-'^im^^) respectively such 
that 

Ai,,{x, = D,{x, + Ceix, + Reix, 0, 

and Ds{x,S,) > 0. We now take the Priedrichs part of the operators Ds{x,Dx) and Ce{x,Dx). Theorem 
A. 3 yields 

Ai^,{x, D^) = DF,e{x, D^) + CfAx, Dx) + Se{x, D^), 

where 5^ has symbol entries in A45-i(r2„), DF,e{x, D^;) is positive and {Ce{x, Dx)u,v) = —{u,Ce{x,Dx)v) 
for all u, V with components in ^(K"). Hence, combining this decomposition with the Sobolev mapping 
properties of generalised pscudodiffcrcntial operators we get the inequality 

(A.3) R^{Ai^s{x,D^)u,u) = {DF,e{x,Dx)u,u) + Re{Ss{x,D^)u,u) > -\\Ss{x,D^)u\\i/2\\u\\_i/2 

> -Ci^e\\u\\-l/2, 

where the net (ci^e)^ depends on the symbol seminorms of the entries of {Ss{x,^))s. The assertion (A.2) 
follows immediately from (A.3). □ 



We conclude this appendix by stating the Garding inequality under different assumptions of moderateness 
at the level of representatives. Focusing on the Fricdrichs part of a pseudodifferential operator one can 
state Theorem A.2 in terms of slow scale symbols and logarithmic slow scale symbols. These two new 
versions, more specific under the point of view of the moderateness assumptions, are stated by replacing 
M with M^'^ and M^^'^ respectively, throughout the statement of the theorem. This leads to the following 
Garding inequalities. 
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Theorem A.5. 



(i) Let A(x, Dx) be an mxm matrix of generalised pseudodijferential operators with symbol in Ggo(^R2ny 
If there exist a representative {Ag)^ of A{x,^) and slow scale nets {c~^)s and {r~^)s such that 

(A.4) ReA,{x,^)>cJ 

for 1^1 > r£ then there exists a slow scale net (ci^e)^ such that 

(A.5) Be{As{x,Dx)u,u) > Cs\\u\\l - ci^e\\u\\-i/2 

for all vectors u with components in =y(]R") and all e € (0, 1]. 

(a) Let A{x, Dx) be an mxm matrix of generalised pseudodifferential operators with symbol in 5go(R2,i) • 
If there exist a representative {As)^ of A{x,^) and logarithmic slow scale nets {cj^)^ and {r~'^)e 
such that 

(A.6) ReA,{x,^)>cJ 
for 1^1 > Te then there exists a logarithmic slow scale net {ci,s)s such that 
(A.7) Re{A,{x,Dx)u,u) > c,\\u\\l - ci,,\\u\W/^ 

for all vectors u with components in ^(M") and all e e (0, 1]. 
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